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PREFACE 


The  Surprisal  Analysis  to  determine  kinetic  information  on  collisional 
processes  developed  by  R.  D.  Levine  and  coworkers  is  the  most  significant 
development  in  molecular  kinetics  in  the  last  decade.  J.  I.  Steinfield  of 
MIT  who  has  implemented  the  surprised  analysis  in  several  molecular  systems 
gave  a series  of  lectures  at  AFWL  while  on  a summer  faculty  term.  This 
inspired  Capt.  C.  Clendening,  J.  I.  Steinfield  and  L.  E.  Wilson  to  apply  the 
Surprisal  Analysis  to  HF  chemical  lasers. 
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SECTION  I 

BACKGROUND  AND  INTRODUCTION 

Detailed  modeling  and  optimization  of  the  HF  (or  DF)  chemical  laser  re- 
quires the  use  of  a large  number  of  rate  constants,  from  several  dozen  to  over 
a hundred  in  the  most  detailed  models.  The  vibrational  deactivation  rates,  in 
particular,  need  to  be  known  over  a wide  range  of  initial  and  final  v states 
and  temperatures.  Only  a handful  of  the  required  rates  are  known  experimentally 
(ref.  1 and  2),  and  it  is  customary  to  guess  those  rates  which  are  not  known. 
Alternatively,  detailed  classical  or  quantum  calculations  of  the  unknown  rates 
can  be  undertaken. 

A more  systematic  approach  to  this  problem  is  now  available  in  the  work 
of  R.  D.  Levine  and  co-workers  (refs.  3 through  22).  These  methods  are  known 
variously  as  the  Thermodynamic  Approach  to  Collision  Processes,  the  Information- 
Theoretical  Analysis,  or  the  Surprisal  Analysis.  The  latter  term,  while  prop- 
erly referring  to  one  specialized  formulation  of  the  general  theory,  is  a con- 
venient term  which  we  shall  use  in  this  report  to  refer  to  this  entire  approach 
and  methodology. 

In  this  report,  we  have  applied  Surprisal  Theory  to  the  synthesis  of  the 
HF-HF  V V',  R,  T deactivation  rate  set.  Detailed  values  for  k(v  -*■  v')  have 
been  generated  for  both  V - (R,  T)  and  V - V deactivation  channels.  This  theory 
also  predicts  the  nature  of  the  temperature  dependence  of  the  HF ( 1 ) + HF(0) 

V - (R,T)  relaxation  rate  between  1 50° K and  4000°K  without  specific  reference 
to  the  details  of  the  intermolecular  potential. 

Among  additional  problems  examined  in  this  report  are  the  use  of  a 
"Populations'  Evolution  in  Time"  (PET)  Approach  to  generate  vibrational  level 
populations  without  explicitly  solving  the  master  equation,  and  Surprisal 
Analysis  of  available  rotational  relaxation  data. 
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SECTION  II 

SURPRISAL  SYNTHESIS  OF  HF-HF  VIBRATIONAL  DEACTIVATION  RATES 
1.  CALCULATION  OF  PRIOR  RATES 

Prior  rates  were  calculated  for  the  Individual  vibrationally  and  rotation- 
ally  inelastic  collisions. 

HF(v,J)  + HF(m,L)  - HF(v',  J')  + HF(m",  L') 

where  v and  m are  vibrational  quantum  numbers  and  J and  L are  rotational  quan- 
tum numbers.  These  have  the  form  (ref.  17) 

k°(vJ,  mL  - vM',  m'L')  = R (2L'  + 1)(2J'  + 1)  eA|A|K1(|A|  )(kT)1/2 


where  Kj(x)  is  the  modified  Bessel  Function  of  the  first  order  and  second  kind, 
and 

. . [<E,  + Em>  * <EJ  + EL>  - <Ev'  * En,'>  + <EJ'  * EL')] 


Since  small  energy  defects  are  important  for  these  calculations,  we  cannot  use 
the  rigid-rotor-harmonic-oscillator  (RRHO)  approximations,  but  instead 


Ev  - (v  + 1/2)  [u)e  - (v  + 1/2)  oigXg] 


Ej  * J(J  + 1 ) [Be  - ae  (v  + 1/2)] 


with  similar  expressions  for  Em  and  EL- 

The  individual  prior  rates  then  need  to  be  averaged  over  rotational  states, 
since  we  wish  to  retain  only  the  vibrational  information.  The  expression  for 
this  is 
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£ £ (2J  + 1)(2L  + l)e-EJ/kTe  ’ EL/kT  k°(vJ'mL  - ^r.n'L') 
k° ( vm  •*  v'm')  - — — — — ■ — - 

( V Ze-Ej/kT  (2J  + l)£e-El/kT  (2L  + 1)  (4) 

J L 

This  k°  would  then  be  multi  pled  by  the  Surprisal  factor,  exp  (-X | A | ) , to  furnish 
values  of  k(vm  -*■  v'm').  In  practice,  X was  found  as  described  below,  and 
k(vm  -*•  v'm')  calculated  as 


k(vtiH-v'm') 


[£e"E0/kT(2J+l)rl  [Se'EL/kT(2L+l)]-1 
J • L 


(2J+1 ) (2L+1  )e”^EJ+EL^kT  k°(vJmL-v'J'm'L')e‘XlEv'+Em'"Ev~Eml 

£ E k°  ( v JmL+v  'J  'm'L ' )s"X  I Ev  '+^m'“Ev_Em I 
J'L'  v'm' 


I 


The  results  of  using  equation  (4)  or  equation  (5)  appear  to  be  numerically 
indistinguishable. 

2.  DETERMINATION  OF  SURPRISAL  PARAMETERS 

The  prior  and  actual  rates  are  related  by  the  "Surprisal"  factor, 

k(vm  -*•  v'm')  = e^"Xo  ” *ilAD  k°(vm  -*■  v'm')  (6) 


If  sufficient  data  are  available,  \l  can  be  found  by  a Surprisal  Analysis 
(refs.  15  through  17);  otherwise,  various  sum  rules  can  be  employed.  Another 
method  relies  on  the  conservation  of  probability.  We  have  calculated  the  rates 
for  the  two  processes  of  "V  - T"  deactivation, 

ki 

HF(1 ) + HF (0 ) ♦ HF(0)  + HF(0) 


and  purely  elastic  scattering  or  exchange, 

k 

HF ( 1 ) + HF (0)  - o HF ( 1 ) + HF(0) 
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or 

HF ( 0)  + HF ( 1 ) 

The  ratio  of  these  rates  depends  only  on  \Ja|.  Since  the  "collision  number" 
for  deactivation  is  known,  namely, 

Zv  • k0/kx  = 75 

we  have  used  this  to  determine  a value  of  (Aj/kT)  = 0.0021  l/(cm_1).  This 
corresponds  to  a At  = 0.438,  or  Ay  = (hccOg/kTjAj  = 8.69.  This  value  appears 
to  be  rather  high  (in  many  systems,  Ay  is  of  order  1),  but  is  consistent  with 
analysis  of  the  HF(v)  + H -*■  HF(v-l)  + H rates,  carried  out  in  reference  17, 
which  gives  = 0.39,  and  with  our  own  analysis  of  (v,0  -*■  v-1 , 1)  rates  in  HF. 

The  absolute  scaling  factor  exp(AQ)  is  found  simply  by  noting  that  the 
measured  value  of  kx (10+00;  300°K)  is  1.0  x 1012  cm3/mole-sec. 

3.  TABULATION  OF  RESULTS 

The  results  of  the  calculation  described  in  parts  (a)  and  (b)  are  74  V-(T,R) 
and  V-V  rate  coefficients  which  appear  to  be  significant  in  the  self-relaxation 
of  HF  levels  v <_  5.  [By  significant,  we  mean  having  a bimolecular  rate  coef- 
ficient _>  1 x 109  cm3/mole-sec.]  These  are  summarized  in  table  1.  The  entire 
calculation  requires  only  10  or  20  minutes  on  a CDC  6600  computer. 

Examination  of  the  individual  rotational  channels  for  the  "V-T"  processes 
f(l)-(9)]  reveals  that  most  of  the  vibrational  exothermicity  is  transferred  to 
nearly  resonant  rotational  levels.  Such  near-resonant  V-R  transfer  was  first 
suggested  by  Houston  et  al.,  (ref.  23),  and  agrees  with  recent  trajectory  cal- 
culations by  Wilkins  (ref.  24).  Multiquantum  (Av>l)  transitions  are  not  sig- 
nificant contributors  to  these  rates,  but  are  important  in  the  V-V  processes 
[(10)- (74) ] in  which  the  total  number  of  vibrational  quanta  is  unchanged  or 
decreases  by  one.  The  individual  V-(R,T)  Av  = -1  rates  do  not  appear  to  scale 
as  v2  or  V2*4,  as  has  been  frequently  suggested;  the  dependence  is  closer  to  v. 
Indeed,  the  basic  sum  rule  [equation  (74),  ref.  17],  namely, 

£ (Ev-  - Ev)k(v  + v';T)  « a[«Ey  »eq  - Ey] 
v' 


8 


AFWL-TR-76-144 


Table  1 

I 

RATE  COEFFICIENTS  FOR  HF-HF  V-(V,R,T)  DEACTIVATION 
PROCESSES,  FROM  SURPRISAL-SYNTHESIS  METHOD 


Process* 

AE  ..** 
vib 

k 

cal/mole 

cm3/mol e-sec 

1. 

HF ( 1 )+HF  - HF(0)+HF 

-11,327 

1.00  x 1012 

2. 

HF(2)+HF  - HF ( 1 )+HF 

-10,834 

2.19  x 1012 

3. 

HF(3)+HF  - HF(2)+HF 

-10,356 

2.88  x 1012 

4. 

HF(4)+HF  - HF(3)+HF 

- 9,889 

3.66  x 1012 

5. 

HF(5)+HF  - HF (4)+HF 

- 9,433 

4.09  x 1012 

6. 

HF(2)+HF  -►  HF(0)+HF 

-22,161 

1.30  x 109 

7. 

HF(3)+HF  - HF ( 1 )+HF 

-21,190 

4.97  x 109 

8. 

HF(4)+HF  -v  HF(2)+HF 

-20,245  . 

9.12  x 109 

9. 

HF(5)+HF  - HF(3)+HF 

-19,322 

1.47  x 1010 

10. 

HF ( 1 )+HF(l ) + HF(2)+HF(0) 

- 493 

4.93  x 1013 

11. 

HF(2)+HF(0)  - HF ( 1 )+HF(l ) 

+ 493 

1.11  x 1013 

12. 

HF(2)+HF(1 ) - HF(3)+HF(0) 

- 971 

3.88  x 1013 

13. 

HF(2)+HF(2)  - HF(4)+HF(0) 

- 1,916 

2.94  x 1013 

14. 

HF(2)+HF(2)  - HF(3)+HF(1) 

- 478 

2.86  x 1013 

15. 

HF(2)+HF(2)  - HF(3)+HF(0) 

-11,805 

4.9  x 1011 

16. 

HF(3)+HF(0)  - HF(2)+HF(1) 

+ 971 

7.97  x 1012 

17. 

HF(3)+HF(0)  •>  HF(1)+HF(1) 

- 9,863 

1.88  x 1012 

18. 

HF(3)+HF(1 ) - HF(4)+HF(0) 

- 1,438 

3.54  x 1013 

19. 

HF (3 )+HF ( 1 ) HF(2)+HF(2) 

+ 478 

5.55  x 1012 

20. 

HF(3)+HF(2)  - HF(5)+HF(0) 

- 2,839 

2.15  x 1013 

21. 

HF (3 )+HF (2 ) - HF(4)+HF(1 ) 

- 945 

2.67  x 1013 

22. 

HF(3)+HF(2)  + HF(4)+HF(0) 

-12,272 

3.22  x 1011 

23. 

HF(3)+HF(3)  - HF(6)+HF(0) 

- 4,116 

1.62  x 1013 

24. 

HF(3)+HF(3)  - HF(5)+HF(1) 

- 1,868 

2.39  x 1013 

25. 

HF(3)+HF(3)  -*■  HF(4)+HF(2) 

- 467 

2.35  x 1013 

26. 

HF(3)+HF(3)  -*■  HF(5)+HF(0) 

-13,195 

1.66  x 1011 

27. 

HF(3)+HF(3)  - HF(4)+HF(1 ) 

-11,301 

5.21  x 1011 

28. 

HF(4)+HF(0)  - HF(3)+HF(1 ) 

+ 1,438 

2.73  x 1012 

29. 

HF(4)+HF(0)  HF(2)+HF(2) 

+ 1,916 

4.48  x 1011 

30. 

HF(4)+HF(0)  HF(2)+HF{1 ) 

- 8,918 

6.14  x 1012 

31. 

HF(4)+HF(0)  - HF ( 1 )+HF(l ) 

-19,572 

6.26  x 109 
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Table  1 (Continued) 

1 * 

RATE  COEFFICIENTS  FOR  HF-HF  V-(V,R,T) 

DEACTIVATION 

PROCESSES,  FROM  SURPRISAL-SYNTHESIS  METHOD 

Process* 

AE,  . ** 

vib 

k 

J 

cal /mole 

cm3 /mol  e-sec 

32.  HF(4)+HF(1 ) - HF(5)+HF(0) 

- 1,894 

3.47  x 1013 

33.  HF(4)+HF(1 ) - HF(3)+HF(2) 

+ 945 

4.12  x 1012 

34.  HF(4)+HF(1)  - HF(2)+HF(2) 

- 9,411 

1.13  x 1012 

35.  HF(4)+HF(2)  - HF(6)+HF(0) 

- 3,743 

1.62  x 1013 

36.  HF(4)+HF(2)  - HF(5)+HF(1) 

- 1,401 

2.69  x 1013 

37.  HF(4)+HF(2)  -*■  HF(3)+HF(3) 

+ 467 

4.32  x 1012 

38.  HF(4)+HF(2)  - HF(5)+HF(0) 

-12,728 

2.42  x 1011 

39.  HF(4)+HF(2)  -»•  HF(3)+HF(1) 

-20,723 

2.3  x 109 

40.  HF (4)+HF (3)  - HF(7)+HF(0) 

- 5,554 

6.72  x 1012 

41.  HF(4)+HF(3)  -►  HF(6)+HF(1) 

- 2,772 

1.92  x 1013 

42.  HF(4)+HF(3)  ♦ HF(5)+HF(2) 

- 923 

2.40  x 1013 

43.  HF(4)+HF(3)  - HF(5)+HF(1) 

-11 ,757 

3.81  x 1011 

44.  HF(4)+HF(4)  + HF(8)+HF(0) 

- 7,335 

2.61  x 1012 

45.  HF(4)+HF(4)  - HF(7)+HF(1) 

- 4,116 

1.20  x 1013 

46.  HF(4)+HF(4)  - HF(6)+HF(2) 

- 1,827 

2.29  x 1013 

47.  HF(4)+HF(4)  - HF(5)+HF(3) 

- 456 

2.25  x 1013 

48.  HF(4)+HF(4)  - HF(6)+HF(1) 

-12,661 

2.12  x 1011 

49.  HF(4)+HF(4)  - HF(5)+HF(2) 

-10,812 

6.54  x 1011 

50.  HF (4)+HF(4)  + HF(5)+HF(1) 

-21 ,646 

1.1  x 109 

51.  HF(4)+HF(4)  - HF(3)+HF(3) 

-19,778 

1.9  x 109 

| . 52.  HF(5)+HF(0)  -*■  HF(4)+HF(1) 

+ 1,894 

8.34  x 1011 

53.  HF(5)+HF(0)  -*  HF(3)+HF(2) 

+ 2,839 

8.43  x 1010 

54.  HF(5)+HF(0)  HF(3)+HF(1) 

- 7,995 

8.68  x 1012 

55.  HF(5)+HF(0)  - HF(2)+HF(2) 

- 7,512 

5.50  x 1012 

56.  HF(5)+HF(0)  -*■  HF(2)+HF(1) 

-18,351 

2.75  x 1010 

57.  HF(5)+HF(1 ) - HF(6)+HF(0) 

- 2,342 

3.28  x 1013 

58.  HF(5)+HF(1)  - HF(4)+HF(2) 

+ 1,401 

1.50  x 1012 

59.  HF ( 5 )+HF ( 1 ) - HF(3)+HF(3) 

+ 1 ,868 

2.51  x 1011 

60.  HF(5)+HF(1)  - HF(3)+HF(2) 

- 8,488 

3.78  x 1012 

61.  HF(5)+HF(2)  - HF(2)+HF(0) 

- 4,631 

1.23  x 1013 

62.  HF(5)+HF(2)  - HF(6)+HF(1) 

- 1,849 

2.80  x 1013 
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Table  1 (Continued) 

RATE  COEFFICIENTS  FOR  HF-HF  V-(V,R,T)  DEACTIVATION 
PROCESSES,  FROM  SURPRISAL-SYNTHESIS  METHOD 

Process*  AEvib** 

cal/mole  cm3/mol e-sec 


63.  HF(5)+HF(2) 

-► 

HF(4)+HF(3) 

+ 923 

3.45  x 1012 

64.  HF(5)+HF(2) 

-> 

HF(3)+HF(3) 

- 8,966 

1.16  x 1012 

65.  HF(5)+HF(3) 

-► 

HF(8)+HF(0) 

- 6,879 

3.55  x 1012 

66.  HF(5)+HF(3) 

-> 

HF(7)+HF(1 ) 

- 3,660 

1.52  x 1013 

67.  HF(5)+HF(3) 

-► 

HF (6 ) +HF ( 2 ) 

- 1,371 

2.53  x 1013 

68.  HF(5)+HF(3) 

-► 

HF(4)+HF(4) 

+ 456 

4.09  x 1012 

69.  HF(5)+HF(4) 

HF(8)+HF(1) 

- 5,441 

6.57  x 1012 

70.  HF(5)+HF(4) 

->■ 

HF(7)+HF (2) 

- 2,715 

1.88  x 1013 

71.  HF(5)+HF(4) 

HF(6)+HF(3) 

- 904 

2.35  x 1013 

72.  HF(5)+HF(5) 

-+■ 

HF(8)+HF(2) 

- 4,040 

1.19  x 1013 

73.  HF(5)+HF(5) 

-► 

HF(7)+HF(3) 

- 1,792 

2.26  x 1013 

74.  HF(5)+HF(5) 

HF(6)+HF(4) 

- 448 

2.22  x 1013 

*In  processes  (l)-(9),  the  vibrational  state  of  the  second  HF 
is  not  explicitly  stated,  although  calculations  indicate  that 
the  rate  of  the  V- (R ,T)  process  does  depend  to  some  extent  on 
HF(v').  Also,  both  forward  and  reverse  processes  are  listed 
separately;  processes  (21)  and  (33),  for  example,  are  of 
course  related  by  detailed  balancing. 

**AEvl-b  is  written  so  that  AEv1-b  < 0 corresponds  to  an  exothermic 

process,  i.e.,  net  conversion  of  vibrational  energy  to  transla- 
tion and  rotation,  and  AEvib  > 0 to  an  endothermic  process. 
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would  seem  to  preclude  vr  scaling  with  r>l . Consider  the  following  simple 
demonstration. 

For  v=l , the  only  significant  rate  is  to  v'  = o,  so  we  have 

(Efl  - EjkO  - 0)  = a[E0  + (E,  - E0)e‘(Ei  ' E°)/kT  - EJ 

- a[(E0  - EJd  - e"hCu)e/kT)] 

The  (Eq-E1 ) cancels,  leaving  the  usual  definition  of  the  vibrational  relaxation 
time, 

_ k(1  + 0) 

a " i-e-hcVkT 

for  v=2,  we  have 

(E0  - E2)  k(2  - 0)  + (£1  - E2)  k(2  - 1)  = a[«Ev  »eq  - E2]  = a(Efl  - E2) 

From  table  1,  we  note  that  k(2-*-0)  is  almost  103  times  smaller  than  k(2+l).  We 
may  thus  write 

(E1  - E2)k(2  - 1)  = a(E0  - E2) 
or 

E - E, 

k(2  - 1)  * F-°-~  F2  k(l  - 0) 

L1  L2 

for  hca)e/kT  » 1.  For  a harmonic  oscillator  this  would  require  k(2-*-l)  = 

2k(l-*-0);  even  including  the  HF  anharmonicity,  we  have  only  k(2-»-l)  = 2.05  k(l-<-0), 
not  4k ( 1-K) ) as  v2  scaling  would  predict. 

We  have  also  noted  that  the  suggestion  that  multiquantum  V-(R,T)  rates 
scale  as  1/Av  (ref.  25)  is  not  borne  out  either.  This  is  compensated,  however, 
by  the  appearance  of  very  rapid  V-V  rates  for  the  higher  v levels. 

The  only  available  experimental  data  with  which  to  compare  these  results  are 
the  flow-tube  experiments  of  Kwok  and  Wilkins  (ref.  26),  which  measure  the  total 
rate  of  deactivation  of  HF(v)  by  HF (0 ) . This  comparison  is  shown  in  table  2. 
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Table  2 

DEACTIVATION  RATES  FOR  HF(v)  by  HF(O) 


k(x  1012 

cm3 /mol  e-sec)* 

V/ 

expt. 
(ref.  26) 

Calculated 

1 

1.0 

1.0** 

2 

12 

13.3 

3 

11 

12.7 

4 

16 

13.0 

5 

35 

19.1 

6 

59 

31 

7 

-- 

45 

*k  is  the  sum  of  V-(T,R)  [e.g.,  HF(3)  + HF(0) 

- HF (2 ) + HF ( 0 ) ] and  V-V  [e.g.,  HF(3)  + HF  (0) 
-►  HF(2)  + HF(1)]  processes. 

**k(l  + 0 -»■  0 + 0)  chosen  for  normalization. 


The  agreement  is  quite  good;  in  particular,  the  slight  decrease  in  overall 
deactivation  rate,  observed  from  (v=2)  to  (v=3),  is  predicted  by  the  Surprisal 
Method.  The  rapid  increase  in  experimental  values  for  v > 3 is  mirrored  by 
the  Surprisal  results,  although  the  rise  is  not  quite  as  fast  as  in  the 
experimental  results. 

Very  recently,  measurements  have  been  carried  out  of  the  v-dependence  of 
the  HF  V-R,  T deactivation  rates  (ref.  27).  These  results  indicate  that  the 
scaling  rule  is  something  like  v1-5,  which  is  much  closer  to  the  Surprisal 
results  than  the  results  of  a v2*4  scaling  would  be. 


14 


. 
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SECTION  III 

TEMPERATURE  DEPENDENCE  OF  k(lCM)O) 


The  temperature  dependence  of  the  deactivation  rate  for 

HF(1)+HF(0)  - HF (0)  HF(0) 

has  been  extensively  studied  (ref.  1).  Its  most  notable  feature  is  an  inverse 
temperature  dependence  below  1000°K,  followed  by  a broad  minimum  which  then 
rises  again  to  an  asymptotic  Landau-Teller  expC-T"1/3)  form  at  high  tempera- 
tures. If  the  Information-Theoretical  analysis  is  to  be  useful,  it  must 
reproduce  this  behavior. 

The  form  of  the  rate  constant  in  equation  (6)  is  derived  from  minimization 
of  the  Lagrangian  (ref.  6). 


+ k(vn»v"nT)  + xV  k(vnH-v"m") 

0 v'm'  vV 


x [(£„.  ♦ E„.)-(Ev  * Em)] 


Detailed  balancing  requires  that  the  reverse  reactions  be  given  by 

k(vm*v'nT)  = e”^Ev'  + Em')  " ^Ev  + Em^kT  k(v'm^vm) 

3 eAE/kT|^v^m^vm)  (9) 

These  reactions  will  become  important  at  high  temperatures,  especially  as  many 
of  the  V-V  processes  have  an  overall  AE  < 1000°K.  Incorporating  equation  (9) 
in  equation  (8)  gives 
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a£~H  k(vm+-v'nT)  log 
^ v'm" 


kjW^m;)  -AE/kTxl 

k°(vm-*v'm')  ' e ' J 


+ - 1)£  k(viTH-v"m")0  + e“AE>/kT) 

0 M'm' 

+ X2  XL  >(vn^v'm')(AE-e‘AE/kTAE) 


v m 


Taking  variations  of  jCwith  respect  to  all  k(vm-*-v"nr)  gives 


(1  + eAE/kT)[log  + U + U0  - 1)(1  + e'AE/kT) 


+ XiAE(l  - e'AE/kT)  = 0 


Dividing  through  by  (1  + e_AE/kT), 

1og  k° (wri+v^m^)  + 1 " 1 + Xo  + Xi  AE  tanh 


Thus,  the  deactivation  rate  constant  is  given  by  the  expressions 


AE 


k(vm-*v'm';T)  = k°(vm-*v'm';T)e'Xo  e"XiAE  tanh  (ZF!^ 


In  general , we  expect  that  both  XQ  and  Xj  inay  be  temperature-dependent 
(ref.  17).  The  former  may  be  eliminated  from  the  calculation  by  using, 
instead  of  equation  (13), 


k(vm.v'n.';T)e-x*4E  t4nh  <AE/2kT> 
k(Vm*V'm'lT)  k0(vnHv'm';T)e-x.At  ta"h 


(10) 


(11) 


(12) 


(13) 


(14) 


Let  us  assume,  provisionally,  that  Xt  is  independent  of  temperature  and  equal 
to  0.0021  l/(cm"1)»  as  found  for  T = 300° K.  We  can  then  calculate  k(10-*00;T) 


from  equation  (14)  for  a range  of  temperatures,  with  AE  = wQ  (HF).  The  results 


of  doing  so  are  shown  in  figure  1.  The  most  striking  observation  is  that  we 
do  find  the  inverse  temperature  dependence  at  low  T,  the  minimum  around  1000°K, 
and  the  rapid  increase  in  k at  high  T.  The  Information-Theoretic  results 


1 
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—0— INFORMATION  - THEORY 
CALCULATION,  EQ(14  ) 

.=^=1  EXPERIMENTAL  VALUES, 
ZZT-)  SEE  FIG.  7,  REF.  1(a) 

RANGE  OF  EXPERIMENTAL  / 
VALUES  AT  300°K  / 


Figure  1.  Temperature  Dependence  of  the  HF(v  = 1)  + HF(v  = 0)  -*•  2HF(v  = 0) 
Deactivation  Rate,  as  Calculated  from  the  Surprisal  Theory.  Also, 
shown  are  several  recent  experimental  measurements  of  this  quantity 
taken  from  reference  1.  The  Surprisal  calculation  is  normalized 
to  the  measured  value  of  the  rate  coefficient  at  300°K. 
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appear  to  have  a deeper  minimum  than  is  found  experimentally,  making  the 
calculated  high- temperature  rates  too  small;  but  this  could  well  be  a result 
of  using  the  for  300°K  over  the  entire  temperature  range.  This  variation 
of  with  T can  be  regarded  as  the  analogue  of  a nonideal  correction  in 
equilibrium  thermodynamics. 

What  is  remarkable  about  this  result  is  that  the  simple  analytic  expression 

in  equation  (13)  or  equation  (14)  makes  no  direct  reference  to  any  details  of 

the  HF-HF  interaction  potential,  such  as  long-range  attractive  forces,  "sticky 
collisions",  hydrogen  bonding,  etc.  It  has  been  customary  to  explain  the 
temperature  dependence  in  figure  1 in  terms  of  such  features.  Of  course,  the 
actual  magnitude  of  the  inelastic  cross  section,  from  which  we  derive  , is 

wholly  dependent  on  the  nature  of  the  potential.  The  point  is,  that  we  have 

reduced  the  description  of  the  problem  to  a single  parameter  which  can  be 
obtained  in  several  different  ways,  and  we  do  not  need  the  complicated  multi- 
parameter  expressions  that  are  often  used  in  calculating  these  rates  (cf. 
ref.  28). 

It  is  also  useful  to  examine  the  high-  and  low- temperature  limiting  forms 
of  equation  (13).  As  T->0,  tanh  (wQ/2kT)-»-l , and  the  rates  are  just  those  which 
would  be  calculated  by  completely  neglecting  the  exothermic  reverse  reactions. 
As  T -*•  °°,  tanh  (u>0/2kT)  T -►  (a)Q/2kT)T  o,  and  ultimately  a purely  statistical 
range,  in  which  k = e"*ok°,  is  attained.  Similar  asymptotic  behavior  has  been 
derived  by  Procaccia  and  Levine  (ref.  17). 


! 
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SECTION  IV 

OTHER  APPLICATIONS  OF  THE  THEORY 


Several  other  aspects  of  Surprisal  metnods,  which  may  be  of  use  in  laser 
kinetics  problems,  have  been  investigated  briefly. 

1.  REMAINING  PROBLEM  AREAS  IN  VIBRATIONAL  SURPRISAL 

Although  the  approach  described  in  the  preceding  sections  has  been  quite 
successful  in  reproducing  the  known  behavior  of  HF-HF  relaxation  rates,  there 
still  remains  a number  of  questions  to  be  explored.  The  most  important  seems 
to  be  the  accurate  formulations  of  sum  rules,  analogous  to  equation  (111)  of 
reference  17  for  V-V'  and  VR-V'R'  exchanges.  This  is  especially  the  case  for 
collisions  between  dissimilar  pairs,  such  as  HF-DF,  HF-H2,  etc.  which  form  an 
important  part  of  real  laser  systems.  In  the  latter  cases,  it  may  be  neces- 
sary to  introduce  as  many  as  four  independent  Surprisal  parameters,  one  for  the 
vibrational  and  rotational  degrees  of  freedom  of  each  species.  Hopefully, 
information  from  one -component  relaxation  experiments  can  be  used  to  determine 
values  of  these  individual  parameters. 

A second  issue  is  that  of  the  temperature  dependence  of  the  X's,  as  dis- 
cussed in  the  previous  section.  In  terms  of  the  thermodynamic  formulation  of 
Surprisal  theory,  treatments  to  date  have  treated  all  reacting  species  as  ideal 
gases.  As  in  equilibrium  thermodynamics,  nonideal  effects  modify  this  behavior. 
Temperature-  and  density-dependent  Surprisal  parameters  may  be  encountered,  and 
it  will  be  necessary  for  the  theory  to  treat  these  effects. 

2.  ROTATIONAL  NONEQUILIBRIUM 

The  importance  of  nonequilibrium  distributions  over  rotational  states  in 
HF  or  DF  laser  active  media  has  recently  been  suggested  by  several  groups 
(refs.  29  through  31).  If  rotational  deactivation  processes  have  to  be 
explicitly  included  In  performance  codes,  the  dimensionality  of  the  problem 
is  measured  to  a frightening  extent.  Instead  of  a few  dozen  rates  connecting 
5 or  6 vibrational  levels,  we  would  have  50  or  60  vibrotational  levels,  with 
some  tens  of  thousands  of  reactions  to  be  considered! 

It  seems  evident  that  Information-Theoretical  methods  will  be  essential 
in  managing  and  coordinating  this  mass  of  data.  A formulation  of  Surprisal 
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Theory  for  rotational  energy  transfer  is  available  (refs.  18  and  19).  We  have 
applied  this  analysis  to  the  available  data  on  rotational  energy  transfer  in 
(refs.  32  and  33).  Hinchen  finds  that  this  data  can  best  be  fit  by  the 
expression  k(ij-*-kl)  «e-ot I ^j“^i  ^+^l~Ek^  with  a = 0.96.  This  is  of  course 
just  the  Surprisal  form  of  the  rate,  k ~ exp (0 | A | /2kT ) , with  0 = 1.92.  While 
more  work  needs  to  be  done  on  that  analysis,  it  does  appear  that  most  of  the 
data  can  be  represented  by  a single  rotational  Surprisal  parameter. 

3.  V-V  CORRELATION  RULES 


Bott  (ref.  34)  has  carried  out  an  empirical  correlation  of  available 
single-quantum  V-V'  energy  exchange,  that  is,  for  processes  of  the  type 


HX(v=l ) + AB(v=0)  - HX(v=0)  + AB(v=l)  + AE 
He  finds  that  the  data  are  best  represented  by 

V "(u,uJ)e-CAE/hv>  (15) 

where  the  energy  transfer  probability  Pyv-  is  defined  as  the  ratio  of  the  trans- 
fer rate  k^  to  the  gas-kinetic  collision  rate;  and  y2  are  the  permanent 
dipole  moments  of  molecules  HX  and  AB;  Vj  is  the  vibration  frequency  of  HX; 
and  AE  is  the  energy  defect. 

It  is  quite  easy  to  see  how  the  main  features  of  equation  (15)  arise  from 
the  Surprisal  form  of  the  rate  constant,  equation  (6).  The  exp(-AE)  dependence 
is  self-evident.  This  is  to  be  multiplied  by  the  Surprisal  factor  Xy.  Now, 
for  hvj  » kT,  we  can  write  X *A  kt/hVj,  where  A is  a constant  of  order  unity 
(ref.  17).  We  thus  obtain  the  exp  (1/hVj)  dependence  as  well.  It  would  be  of 
interest  to  carry  out  a more  rigorous  derivation  of  this  result,  particularly  to 

y 

see  how  the  (uj)j2)  dependence  enters  into  the  pre-exponential  (Ay  or  e o) 
factor. 


4.  POPULATIONS'  EVOLUTION  IN  TIME  APPROACH 

Recently,  Procaccla  «nd  Levine  (ref.  20)  have  indicated  a method  whereby 
the  populations  In  a multi-level  system  relaxing  toward  equilibrium  may  be 
explicitly  calculated  without  reference  to  detailed  rate  constants.  In  this 
"Populations'  Evolution  in  Time  (PET)"  Approach,  the  populations  are  given  as 
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nv(t)  = e"e^Ev  " \><*>  { 

The  time- dependent  parameter  @ is  found  by  using  a type  of  sum  rule,  which 
requires  that 

R[0(t)]  = <V(t)>  I 


<V(t)  > = < V(~)  > + [<  V(tQ)  > - < V (“>)  >]e"^t_to^T 


EE,e 

. V 


-6  ( t ) E 


The  normalization  parameter  XQ(t)  is  then  found  simply  as 


Kit)  = In  (2>‘e(t)  Ev 


Thus,  the  only  inputs  required  are  the  initial  departure  of  the  system  from 
equilibrium  at  time  tQ,  [<  V(tQ)  > - < V(«)  >],  and  the  bulk  relaxation  time  t. 

We  have  applied  the  PET  method  to  a simple  model  of  relaxation  in  pure  HF 
following  an  injection  of  population  at  some  specified  v",  such  as  would  occur 
in  pumping  av'  + O transition  with  a pulsed  laser.  A full  master  equation 
(for  0 _<  v ^ 5),  using  the  deactivation  rate  constants  in  table  1,  was  solved 
by  numerical  integration.  Equations  (16)  through  (19)  were  then  solved  for 
the  same  initial  conditions,  and  with  a relaxation  time  t consistent  with  the 
rates  used  in  the  model.  The  comparison  is  shown  in  figure  2.  We  see  that  the 
PET  method  gives  essentially  perfect  agreement  with  the  full  master  equation 
for  each  vibrational  level,  as  had  been  previously  noted  (ref.  20).  The  only 
quantity  that  is  not  well  determined  is  the  relative  proportions  in  higher 
v states  following  the  excitation  pulse,  which  is  determined  by  V«-**V  exchanges. 
Levine  suggests  that  these  can  also  be  calculated  by  adding  an  additional  VV 
constraint  to  equation  (16)  to  give 
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nv(t) 


re(t)Ev+X(t)|Ev  - E.|-X0(t) 


where  v=i  is  the  level  initially  displaced  from  equilibrium. 

On  a more  profound  level,  major  simplifications  may  be  possible  in  the  way 
in  which  we  model  complex  chemical  laser  systems,  beyond  merely  generating  sets 
of  rate  constants.  When  using  the  modeling  codes,  after  all,  we  are  not  really 
interested  in  the  detailed  values  of  these  rate  constants;  the  information  we 
actually  require  is  [n(x,t)],  the  time  and  space  variation  of  the  concentrations 
of  particular  species,  such  as  HF(v).  We  need  the  rate  constants  only  as  input 
to  a kinetic  master  equation  which  must  then  be  coupled  to  optics  and  gas- 
dynamic  equations  to  find  these  concentrations.  The  whole  idea  of  the  Infor- 
mation-Theoretic approach  is  to  find  expressions  only  for  the  data  we  want  to 
know,  while  avoiding  detailed  consideration  of  large  masses  of  unnecessary  data. 
It  may  be  possible  to  apply  the  PET  approach  to  the  gas-dynamic  chemical  laser 
problem,  in  an  attempt  to  find  a direct  route  to  [n (x ,t ) ] while  bypassing  as 
much  of  the  detailed  kinetics  as  possible.  This  problem  has  been  treated  in  a 
general  way  by  Ben-Shaul  and  Hofacker  (ref.  22),  and  further  analysis  of  specific 
systems  weuld  be  highly  desirable.  If  it  were  indeed  possible  to  streamline 
performance  codes  by  eliminating  much  of  the  detailed  kinetics,  a much  more 
thorough  modeling  of  the  optics  and  gas-dynamic  parameters  would  then  be  both 
possible  and  practical. 


AFWL-TR-76-144 


t 

( 


i 


i 


SECTION  V 
CONCLUSIONS 

Our  initial  experience  in  applying  Surprisal  Theory  to  the  HF  relaxation 
problem  has  been  that  a great  deal  of  insight  may  be  obtained  with  the  theory 
as  it  exists,  but  that  a number  of  questions  remain  to  be  answered.  These 
questions  pertain  both  to  fundamental  aspects  of  the  theory,  its  range  of 
validity,  etc.,  and  to  practical  aspects  of  its  implementation.  Among  the 
problems  that  need  to  be  investigated  are  the  best  methods  for  finding  the 
\ parameter,  proper  formulation  of  sum  rules  for  coupled  vibration-rotation 
systems,  and  so  on.  Asymptotic  forms  of  the  rate  constants  for  limiting  hi gh- 
and  low- temperature  t would  also  be  useful. 

The  Surprisal  Theory  has  at  this  point  received  an  impressive  number  of 
validations,  both  from  experimental  and  computed  data  bases,  in  HF  and  many 
other  systems.  If  we  accept  the  point  of  view  that  this  theory  is  an  essen- 
tially correct  description  of  Nature,  then  our  entire  approach  to  state-to-state 
chemical  kinetics  needs  to  be  revised.  The  present  approach  is  to  carry  out 
ever  more  sophisticated  and  detailed  calculations,  and  to  validate  these  with 
ever  more  marginal  experiments.  The  Information-Theoretical  approach  states 
that  essentially  all  microscopic  rate  data  for  a given  system  is  interrelated, 
and  determined  by  one  or  two  parameters.  The  purpose  of  an  experiment  is  then 
to  establish  the  value  of  these  parameters  by  very  careful  measurements  of  a 
basic  system  property,  such  as  the  relaxation  time.  A calculation  is  deemed 
to  be  correct  once  it  is  done  carefully  enough  to  reproduce  the  Surprisal- 
synthesized  results,  and  all  theoretical  methods  — classical  trajectories, 
close-coupled  quantum  states,  R-matrix  theory,  etc.  — should  converge  to  these 
results  if  properly  done.  In  applying  kinetic  information  to  the  analysis  and 
modeling  of  real  systems,  Surprisal  Methods  again  offer  the  prospect  of  much 
more  rapid  and  efficient  modeling  procedures  than  are  presently  being  used,  and 
also  make  it  possible  to  bypass  much  of  the  detailed  information  now  required. 
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